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Abstract 

Cronobacter spp. {Enterobacter sakazakii) is an important pathogen contaminating powdered infant formula (PIF). To 
describe the genotypic diversity of Cronobacter isolated in China, we identified the isolates using fusA allele sequencing, and 
subtyped all of the isolates using pulsed-field gel electrophoresis (PFGE), multi-locus sequence typing (MLST), and multiple- 
locus variable-number tandem-repeat analysis (MLVA). A total of 1 05 isolates were identified, which included C. sakazakii (58 
isolates), C. malonaticus (30 isolates), C. dublinensis (1 1 isolates), C. turicensis (5 isolates), and C. muytjensii (1 isolate). These 
isolates were showed to have 85 PFGE-patterns, 71 sequence types (STs), and 55 MLVA-patterns. Comparisons among the 
three molecular subtyping methods revealed that the PFGE method was the most distinguishable tool in identifying clusters 
of Cronobacter spp. through DNA fingerprinting, and MLST method came second. However, ESTR-1, ESTR-2, ESTR-3, and 
ESTR-4 were not effective loci for subtyping Cronobacter spp. such that the MLVA method requires further improvement. 
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Introduction 

Infants fed with powdered infant formula (PIF) contaminated 
with Cronobacter (Enterobacter sakazakii) can suffer from necro- 
tising enterocolitis, bacteraemia, and neonatal meningitis [1—4]. 
Although these infections are relatively rare overall, they are of 
considerable significance, because infants comprise the highest-risk 
age group. While these infections are lethal for infants, recent 
studies reported possible diarrhea, wound infections, and urinary 
tract infections among immunocompromised people, particularly 
the elderly [5-8] . While the reservoir of Cronobacter is unknown in 
many cases, PIF was reported as a key source for infantile 
infection. Since 2004, the Food and Agriculture Organization of 
the United Nations (FAO) and the World Health Organization 
(WHO) have addressed the issue of E. sakazakii in PIF: 
accordingly, many studies of the phenotype and genotype of 
Cronobacter were conducted. In particular, Cronobacter was 
proposed as a new genus for E. sakazakii, and the genus consists 
of ten genomospecies [9-11]. These bacteria are distributed 
throughout nature and are found in plants, spice, herbs, water, dry 
foods, and starches [12-15]. Although there were no clinical 
infections reported, these strains have been widely isolated from 
food and from the environment in China. 

Several molecular techniques have been developed for subtyp- 
ing Enterobacter sakazakii, such as ribotyping, random amplifica- 
tion of polymorphic DNA, pulsed-field gel electrophoresis (PFGE), 
multiple-locus variable-number tandem-repeat analysis (MLVA), 
and multi-locus sequence typing (MLST) [16-18]. PulseNet 



International, which is the international molecular subtyping 
network for foodborne disease surveillance, considers PFGE as a 
gold-standard method for bacterial molecular typing. MLVA is the 
second generation of molecular subtyping method. MLST data is 
comparable across laboratories worldwide via the central curated 
database(http://pubmlst.org/cronobacter). To compare the geno- 
types of Cronobacter isolates from different sources, years and 
regions in China, we identified 105 Cronobacter isolates obtained 
from 2006 to 2011 from 11 Chinese provinces, and using PFGE, 
MLVA and MLST method to type these isolates. We also 
investigated which molecular subtyping method is more effective 
in determining the genetic relationship of Cronobacter isolates for 
the purpose of surveillance and outbreak investigation. 

Materials and Methods 

Bacterial strains 

From 2006 to 2011,105 isolates were collected from 11 
provinces in China, and their sources included infant food, 
potable water, and rectal swabs of healthy humans. Three rectal 
swabs from healthy humans were collected by Jinan municipal 
CDC, and the participants were informed and had provided their 
written consent to participate in this study, which was approved by 
the Jinan municipal CDC Ethics committee. As controls, the 
isolates ATCC 51329 and ATCC 29544 were provided by the 
Chinese Academy of Inspection and Quarantine (CAIQ). All 
isolates were confirmed by biochemical test using kit API 20E 
(bioMerieux, Marcy l'Etoile, France) and real-time PCR using a 
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primer set and probe targeting the dnaG gene on the macromo- 
lecular synthesis operon [19]. Isolates were cultured on TSA agar 
(Oxoid, Basingstoke, United Kingdom), and they were cryopres- 
ervation at — 80°C for long-term storage. 

PFGE 

Cronobacter isolates were typed as described in the PulseNet 
standardised PFGE Cronobacter subtyping protocol [1 7] . Xbal and 
Spel were the primary and secondary restriction enzymes, 
respectively. DNA fragments were separated by electrophoresis 
(CHEF Mapper, Bio-Rad Laboratories, Hercules, California, US) 
through a 1 % (w/v) agarose gel (Seakem Gold, Rockland, Maine, 
US) in 0.5 x TBE buffer at 6 volts/ cm using an initial switch time 
of 1.8 s and a final switch time of 25 s. Gels were stained in 
deionised water containing GelRed Nucleic Acid Gel Stain 
(Biotium, CA, US), and visualised under U.V. light using a 
GelDoc XR + system (Bio-Rad laboratories, Hercules, California, 
US). Xfel-digested Salmonella serotype Braenderup H9812 was 
used as the molecular weight standard. Dendrograms were created 
using BioNumerics software version 5.1 (Applied Maths, Sint- 
Martens-Latem, Belgium) using the DICE coefficient, the 
unweighted pair group method with arithmetic means (UPGMA), 
and a band position tolerance of 1.5%. 

Genomic DNA Isolation 

Bacteria were cultured on TSA agar at 37°C for 18-24 h, and 
an isolated colony was inoculated into 5 ml tryptone soy broth 
(Oxoid, Basingstoke, UK). Bacteria from 1 ml of overnight culture 
was recovered by centrifugation at 1 0,000 xg for 2 min. Total 
DNA was prepared using the QIA amp DNA Mini kit (Qiagen, 
Maryland, US) and quantified using a spectrophotometer (IM- 
PLEN, Gemany) for MLST and MLVA method. 

MLST 

The seven multilocus loci of Cronobacter were atpD,fusA, glnS, 
gltB, gyrB, infB, and ppsA. The primers and the detailed 
amplification conditions for these fragments were available at 
http://pubmlst.org/cronobacter/info/protocol.shtml; Baldwin et 
al. 2009. The amplified fragments were purified using a PCR 
purification Kit (Qiagen, Maryland, US) according to the 
manufacture's instructions. DNA sequencing was performed with 
an ABI 3730XL DNA Analyzer (Applied Biosystems, CA, US). In 
order to identify Cronobacter species, the phylogenetic and 
molecular evolutionary analyses of fusA allele were conducted 
by using MEGA (Molecular Evolutionary Genetics Analysis) 
version 5.0 [20]. Concatenated DNA sequences for sequence 
types not detected in the present study were downloaded from the 
MLST database. 

MLVA 

Four loci, ESTR-1, ESTR-2, ESTR-3, ESTR-4 [16] were 
amplified using the primers at a concentration of 0.2 uM with EX 
Taq DNA polymerase (TBI, Dalian, China) in a 20 (0.1 final 
reaction volume using the manufacturer's reaction buffer. 
Following amplification, all PCR products were electrophoresed 
using an Applied Biosystems 3730XL DNA Analyser (Applied 
Biosystems, CA, US). The band intensity of each amplicon was 
examined prior to analysis with Gene Mapper 4.0 software 
(Applied Biosystems, CA, US). All four loci were amplified three 
times from each isolate to ensure reproducibility. To confirm the 
numbers of alleles at the four loci, PCR products were processed 
using the BigDye kit and sequenced on an ABI 3730XL capillary 
DNA analyser (Applied Biosystems, CA, US). 



Results 

Species identification 

Based on Joseph's study [18], the 105 isolates in this study were 
identified as follows: C. sakazakii (58 isolates), C. malonaticus (30 
isolates), C. dublinensis (11 isolates), C. turicensis (5 isolates), and 
C. muytjensii (1 isolate) (Fig. 1). In addtion, the complete 16S 
rRNA gene was also sequenced and additional biochemical tests 
were performed, including the indole test, acid production from 
dulcitol and methyl-a-D-glucopyranoside, and malonate utiliza- 
tion (data not showed). 

PFGE 

BioNumerics software analysis showed the 105 isolates demon- 
strated 85 distinguishable Xbal-Spel PFGE patterns, and the 
control strains (C. sakazakii ATCC29544, C. muytjensii 
ATCC51329) showed the two PFGE patterns different from the 
isolates in this study (Fig. 2). A high degree of genetic diversity was 
revealed using PFGE after genomic DNA digestion, and the 
discriminatory index was 0.9940, as it was calculated using 
Simpson's diversity index. Nevertheless, when the isolates digested 
with Xbal and Spel respectively, 83 and 86 PFGE patterns were 
demonstrated. So the identification rate of PFGE with Spel may 
be higher than with Xbal. 

MLST 

The seven defined alleles of each isolate were sequenced, and 
each house keeping gene atpD, fusA, glnS, gltB, gyrB, infB, and 
ppsA was found to possess 14, 5, 15, 25, 24, 14, and 27 new types 
of allele sequences, respectively. All the sequences were submitted 
to the MLST database and assigned new numbers, which were 
showed in Fig. 1. The 105 strains were in 71 STs, of which only 15 
types had been reported by other countries. Although many new 
STs were found, these STs belonged to previous clonal complex 
(CC, single or double locus variant of dominant sequence type), 
such as ST183, 203, 235, 241, 242 belonged to ST4; and ST159, 
161, 201, 238 belonged to ST7; and so on. The isolate C. sakazakii 
ATCC29544 showed the same ST with one isolate (BJCR0014), 
whereas C. muytjensii ATCC51329 showed different ST from all 
the isolates in this study. Collectively, the above results suggest 
Cronobacter isolated from China had a high diversity, particularly 
for the isolates identified to C. dublinensis and C. turicensis. 

MLVA 

Using the MLVA method developed by Mullane el al., we 
calculated the numbers of alleles for the ESTR-1, ESTR-2, 
ESTR-3, and ESTR-4 loci: they were 6, 4, 12, and 2, respectively. 
For certain isolates, there was no PCR amplicon at a variable 
number of tandem repeats (VNTR) locus, which was labeled as 
"NP" meaning"no product", and the amplicon was absent in 1 1% 
of the isolates at the ESTR-1 locus, 17.1% at the ESTR-2 locus, 
49% at the ESTR-3 locus, and 22% at the ESTR-4 locus. The 
genetic diversity of VNTRs for each locus was calculated using 
Simpson's diversity index, and the values for ESTR-1, ESTR-2, 
ESTR-3, and ESTR-4 loci were 71.68%, 61.16%, 89.55%, and 
32.67%, respectively. Although the polymorphism indicator at the 
ESTR-3 locus was the best, the rate of missing amplicons was very 
high. For VNTRs in ESTR-1 locus, the upstream and down- 
stream sequences were not conserved in 12% of the isolates. To 
distinguish these results from conserved sequences, the repeat copy 
number was determined by negative correlation (Fig. 2). Minimum 
spanning tree (MST) revealed that the 105 isolates were divided 
into 55 MLVA patterns, and the control strains (ATCC29544, 
ATCC51329) showed two MLVA patterns different from the 
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Figure 1. Maximum likelihood tree of the fusA alleles and MLST profiles for the differentiation of Cronobacter species in this study. 

doi:10.1371/journal.pone.0102179.g001 
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Figure 2. Dendrogram combining PFGE patterns of Xba\- and Spel- digested DNA from Cronobacter spp. in this study. Isolates 
information, STs, CC and MLVA allele string are also included. 
doi:1 0.1 371 /journal.pone.01 021 79.g002 
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isolates in this study. The isolates belonging to the different species 
were gathered in one cluster (Fig. 3) 

Discussion 

Iversen et al. used f-AFLP fingerprints, ribopatterns, full-length 
16S rRNA gene sequences, and biochemical tests for species 
typing in previous study [9]. While 16S rRNA gene sequences 
failed to distinguish C. sakazakii from C. malonaticus, and lack of 
continuity between Cronobacter biotypes and species was reported 
[21]. Here the isolates' species identification was determined by 
fusA allele sequence typing based on Joseph's study. Similar to 
other studies [22,23], C. sakazakii was isolated from all sources, 
and was the most common Cronobacter species. Except for four 
isolates from food sources, the other 12 isolates, which were 
classified as C. dublinensis or C. turicensis, were from environ- 
ment. This indicats that different species of Cronobacter is 
prevalent in different econiche. In addition, of the three isolates 
from rectal swabs of healthy individuals, two were identified as C. 
sakazakii (ST40, CC45; ST 154), and one was identified as C. 
malonaticus (ST161, CC7). It has been reported that the 
Cronobacter spp. that cause disease were mainly C. sakazaikii 
and C. malonaticus. These results support the hypothesis that C. 
sakazakii and C. malonaticus may be associated with humans 
while other species with plants and the environment [24] . 

PFGE is well-established and widely used for the molecular 
subtyping of bacteria, including Enterobacter sakazakii [25-28]. 
Brengi et al. modified the PulseNet standard protocol in 2012, and 
we performed and evaluated its applicability for analyzing the 
isolates from China. Most isolates with the same Xbal-Spel PFGE 
pattern were isolated from the same source, region and year, 
except for one cluster (PFGE-pattern X71S33). The strains of this 
cluster were from 5 provinces, and isolated from 2005 to 2008, 
which maybe demonstrate this type of isolates was prevalent at 
that time in China. We also found that the cluster (PFGE-pattern 
X67S26) was existent in Sichuan Province from 2010 to 2011. 
Unexpectedly, the fingerprint similarity rate from the same species 
was not high. Although several isolates demonstrated the same 
PFGE pattern when they were digested with Xbal, they showed 
different PFGE patterns when using Spel, and the converse was 



true. Therefore, it was important to perform PFGE on Cronbacter 
using at least two restriction enzymes as a useful tool for 
surveillance and outbreak investigation. 

The fusA gene was more conservative comparing with other 
house-keeping genes, only five new alleles were found in 105 
isolates, so it was more suitable for species. While, other genes 
possessed more than ten types of new alleles, with the result that 
many novel STs, which were quite different from the alleles 
reported by other countries, were found in this study. In the same 
STs, some parts of isolates from the same region and time, such as 
ST155, ST17, ST201, ST 240, and ST247; however, the isolates 
of ST7, ST4, ST1, ST23, ST184, and ST188 were from different 
source, region and time. MLST method was not enough 
distinguishable for subtyping Cronobacter spp., but it was a 
convenient method to compare the data across laboratories 
worldwide. For example, C. sakazakii CC4 was reported to be 
associated with neonate meningitis by other countries [29], while it 
was found in food and water in China. There were 3 isolates from 
healthy individuals in China were identified to C. sakazakii ST40 
(CC45), C. sakazakii ST 154, and C. malonaticus ST161(CC7), 
which have not been reported to associated with clinical cases. 

MLVA-Pattern 28 had the most number of isolates, while the 
isolates from this cluster were from the different regions, sources, 
and year. The identification rate of MLVA method was not high. 
Mullane et al. obtained similar results in their study, and they 
attributed this to the C. sakazakii strain ATCC BAA-894 VNTR- 
based primers, which were designed prior to the Cronobacter 
genus definition. In our study, it was found that the isolates that 
failed to amplify did not always belong to C. sakazakii. Therefore, 
this MLVA method was limited to the full genus as currendy 
recognized, and should be further optimized for effective 
subtyping of Cronobacter. 

Comparing the three molecular subtyping methods by the 
Simpson's diversity indexes, the D-value of PFGE method was the 
highest for each species (Table 1). Although certain isolates with 
the same PFGE-pattern were not in same STs, they were in same 
CCs, while the opposite results were not obtained; certain isolates 
with similar PFGE-Patterns were not gathered in one cluster using 
MLVA methods: the opposite results also were obtained. 




Figure 3. WIST of MLVA analysis of 107 isolates of Cronobacter spp. in this study. Different species were labelled with different colour: C. 
sakazakii with green, C. malonaticus with red, C. dublinensis with purple, C. turicensis with yellow, and C. muytjensii with blue. 
doi:1 0.1 371 /journal.pone.01 021 79.g003 
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Table 1. Numbers of patterns and D-values of PFGE, MLST and MLVA among different Cronobacter species. 







PFGE 




MLST 




MLVA 




No. of patterns 


D 


No. of patterns 
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No. of patterns 
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C. sakazakii (n = 58) 


52 


0.9958 


41 


0.9806 


40 


0.9831 


C malonaticus (n = 30) 


16 


0.8782 


13 


0.8207 


8 
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C. dublinensis (n = 1 1 ) 


11 


1 .0000 


11 


1 .0000 


8 


0.8909 


C turicensis (n = 5) 


5 


1.0000 


5 


1 .0000 


1 


0.0000 


C. muytjensii (n = 1 ) 



doi:1 0.1 371 /journal.pone.01 021 79.t001 



Conclusions 

In conclusion, there were polymorphisms between the Crono- 
bacter spp. isolated in China. The PFGE method was a useful tool 
in identifying clusters of Cronobacter spp. through DNA finger- 
printing, and assisted in retrospective investigations to identify 
food and environmental sources of Cronobacter . Although the D- 
value of the MLST method was lower than the PFGE method for 
each species, it was a convenient method to compare the data 
across laboratories worldwide. However, ESTR- 1 , ESTR-2, 



ESTR-3, and ESTR-4 were not effective loci for subtyping the 
Cronobacter spp. from China. And new loci are required to 
establish an MLVA-based protocol in the future. 
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